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ABSTRACT 

Aims. W49A is a giant molecular cloud which harbors some of the most luminous embedded clusters in the Galaxy. However, the 
explanation for this starburst-like phenomenon is still under debate. 

Methods. We investigated large-scale Spitzer mid-infrared images together with a Galatic Ring Survey 13 CO 7=1-0 image, 
complemented with higher resolution (~ 1 1") 13 CO 7 = 2-1 and C ls O 7 = 2-1 images over a ~ 15x13 pc 2 field obtained with the 
IRAM 30m telescope. 

Results. Two expanding shells have been identified in the mid-infrared images, and confirmed in the position-velocity diagrams made 
from the 13 CO 7 = 2-1 and C ls O 7 = 2-1 data. The mass of the averaged expanding shell, which has an inner radius of x 3.3 pc and 
a thickness of as 0.41 pc, is about 1.9 x 10 4 M©. The total kinetic energy of the expanding shells is estimated to be ~ 10 45 erg which is 
probably provided by a few massive stars, whose radiation pressure and/or strong stellar winds drive the shells. The expanding shells 
are likely to have a common origin close to the two ultracompact Hn regions (source O and source N), and their expansion speed is 
estimated to be ~ 5 km s _1 , resulting in an age of ~ 3-7 x 10 5 years. In addition, on larger (~ 35 x 50 pc 2 ) scales, remnants of two gas 
ejections have been identified in the 13 CO 7=1-0 data. Both ejections seem to have the same center as the expanding shells with a 
total energy of a few times 10 50 erg. The main driving mechanism for the gas ejections is unclear, but likely related to the mechanism 
which triggers the starburst in W49A. 
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1. Introduction 

A variety of processes have been suggested that might trigger the 
formation of stars, namely large-scale supernova explosions or 
small er-scale stellar winds, outflow shocks, or io nization fronts 
(see iMcKee & Ostrikerl 120071 lElmegreerl Il998). Possible trig- 
gers are newly formed massive stars, which provide energetic 
feedback to the interstellar medium (ISM) in the form of ex- 
panding Hn regions or strong stellar winds. A recent study of 
the ubiquitous dust bubbles observed in the mid-infrared (MIR, 
5^-0 fim) in the Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE; Benj amin etafl 120031) shows that, 
given the positional coincidence with known Hn regions, many 
MIR bubbles are indeed produced by hot young stars, and are 
likely driven by exp anding Hn shells (IChurchwell et all 2006; 
■ IWatson et al.. 2008). In addition, for reasonable assumptions on 
the ambient density, t he numerical simulations of expandi ng Hn 
regions carried out bv lHosokawa & Iriu tsuka (2005, 2006) agree 
with some observations, e.g.. lWatson et al.1 (I2008I) . 

The giant molecular cloud (GMC) associated with the pow- 
erful radio continuum source W49A is part of the W49 complex 
which also inc ludes the supernova remnant W 49B. At a distance 
of 11.4 kpc dGwinn. Moran. & Reidl |1992|) . W49A has an in- 
frared lumino s ity of > 10 7 L0 and a total mass of ~ 10 6 Mq 
dSievers et all Il99ll) . and contains over a dozen of ult r acom- 
pact Hn (UC Hn) r egions dDe Pree. Mehringer. & Gossl 119971; 
iDe Pree et all l2000h as well as the brightest water maser cluster 
in ou r Galaxy (~1 Lq alone in the 22 GHz line. iGenzel et all 
fl978h . The embedded massive stellar cluster found in W49A 
indicates a somewhat more evolved stage of massive star for- 



mation in this region ves & Homeierl l2003h . Many authors 
have tried to explain why W49A, a Galactic mini-starburst, has 
forme d so many massive stars in such a short time. IWelch et"aL1 
(1 19871) concluded the existence of a large-scale inside-out gravi- 
tational collapse based on their interpretation of the HCO + spec- 
tra measured to ward the ring-like configuration of UC Hn re- 
gions. However, Serabv nT Glisten. & Schulzl d 19931) argued that 
the two components of the double-peak line profile seen in CS 
and C 34 S come from different clouds, and suggested that mas- 
sive star formation in W49A is triggered by a large-scale cloud- 
cloud collision. Recently, it has been argued that one cannot 
distinguish between a global collapse model of a large cloud 
and a multiple-cloud model, given the existing HCO + and C ls O 
data (Williams. Dickel. & Auerl 12004). In this paper, we present 
new evidence of expanding shells in W49A delineated by dense 
clumps within these shells. The expanding shells provide a nat- 
ural explanation for the observed molecular line profiles, and 
triggered massive star formation in W49A. 



2. Observations 

Observations of the 13 CO J = 2 - 1 and C I8 7 = 2-1 tran- 
sitions toward W49A were carried out in 2006 November and 
2007 November using the IRAM 30 m telescope on Pico Veleta 
in Spain. The focus was checked on Uranus, and the pointing 
was checked on the strong continuum emission from the BL Lac 
object 1749+096 and Uranus every hour. The accuracy of the 
pointing was ~ 2"— 3" with a telescope tracking accuracy of 1". 
The chopper wheel method was used to calibrate the data. 
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Fig. 1. The Spitzer-GLIMPSE IRAC 4.5 /mi and 8.0 //m images of W49A. Black triangles denote the UC Hn regions observed by 
iDe Pree.M ehringer, & Gossl d!997l) . The black dashed ellipses are the suggested expanding shells (shells 1 and 2), and the black 
filled star represents the center of the expanding shells and ga s ejections (see text). The larger black unfilled star represents the 
embedded cluster center according to lAlves & Homeierl (120031) . The white circle in the left image indicates the maximum radius 
(5.3 pc) used in the circular PV diagram shown in Figure[3] and the two white dashed lines represent the linear PV cuts (Fig. [3}. 
Red lines and arrows in the right image indicate the directions of gas ejections. 



We used the H Eterodyne Receiver Array, HERA 
dSchuster et al.L 120041) . to obtain large On-The-Fly (OTF) 
maps of W49A. HERA is a nine-pixel receiver with dual 
polarizations in a square center-filled 3x3 array, and can be 
tuned to different sky frequencies with different backends 
attached. In particular, simultaneous observations of 13 CO 
7 = 2-1 and C ls O 7 = 2-1 are possible. A dewar rotation 
angle of 9.6°of HERA was chosen to have Nyquist-sampled 
maps. Maps sizes of ~ 4f5x4!0 were obtained for the 13 CO and 
C ls O 7 = 2-1 lines with an HPBW of ~ 11". Every scan in 
right ascension was followed by one in declination in order to 
reduce striping patterns. The calibration was done for each OTF 
subscan and reduced in the MIRA software of the GILDASQ 
package. OFF scans on the reference position (-600", -600") 
were done before and after each OTF subscan. The Versatile 
SPectrometer Array (VESPA) was chosen as the backend for 
HERA configured to a bandwidth of 160 MHz for each receiver 
element and 0.32 km s _I resolution. The data were reduced 
using the standard procedures in the GILDAS package, and 
were corrected to main beam temperature Tmb units from the 
antenna temperatures T* A (i.e., Tmb = T* a It]mb) which were only 
corrected for atmospheric absorption, rear spillover, and ohmic 
losses. The efficiency, 77MB. is given by the ratio between the 
main beam efficiency (named B e ff; it is 0.51 at 220 GHz in the 
IRAM documentation^)), and the forward efficiency F e g (0.91 at 
220 GHz). 



1 http://www.iram.fr/IRAMFR/GILDAS 

2 http://www.iram.fr/IRAMES 



3. Results 

Figure Q] shows evidence of shell structures in W49A in the 
4.5 /mi and 8.0 /mi images from the GLIMPSE archive. They 
were obtained wi th the Spitzer Infrared Array Camera (IRAC; 
iFazio et all 120041) . The 4.5 /mi band, which contains a pure ro- 
tational transition of H2, mainly traces the distribution of warm 
molecular hydrogen. The 8.0 //m band, which includes two aro- 
matic infrared features at 7.7 /mi and 8.6 /mi, is believed to be 
dominated by the emission of polycyclic aromatic hydrocarbon 
(PAH) molecules primarily excited by UV radiation. The angu- 
lar resolution of IRAC ranges from ~ l'/5 at 3.6 //m to ~ Y.'9 
at 8.0 /mi. In two other IRAC bands (3.6 /mi and 5.8 /im), how- 
ever, the same structures can also be clearly seen. The shells 
were first identified by eye and then investigated in the molecu- 
lar data cube. Two MIR shells have been identified (Fig.[TJ with 
long axes of 55"-58" and short axes of 28"-35", which implies 
an average size of ~ 2-3 pc in radius. The average thickness of 
the MIR shells is estimated to be 7" 5 + 1'.'9 (0.4±0. 1 pc) from fits 
of Gaussian profiles to the shells. These two MIR shells lie in the 
east-west and north-south directions, and have a similar size and 
elliptical shape. In addition, they seem to share the same center. 

However, it is difficult to identify any similar shell structure 
in the channel maps of the 13 CO and C ls O 7 = 2-1 lines (see 
Fig. |2| because of the lower angular resolution and the com- 
plexity/dumpiness of the molecular emission. Position-velocity 
(PV) diagrams are thus applied to investigate this shell structure, 
and are used to examine whether a shell appears as an expand- 
ing ring/bubble in the molecular line data. Linear PV diagrams 
of 13 CO 7 = 2 - 1 and C ls O 7 = 2 - 1 are shown in Figure^ cor- 
responding to the east-west cut (Cut I) and the north-south cut 
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Fig. 2. The IRAM 30m HERA channel maps of 13 CO 7 = 2-1 (images) overlaid with the C 18 7 = 2-1 emission in black contours 
with a resolution of 1 1". The contours are plotted from 1 to 15 K in steps of 2 K. The white ellipses represent two suggested shells, 
and the white star marks the shell expansion center. The (0", 0") position is located at the region around source G. 



(Cut II) in Figure Q] where the 13 CO and C 18 counterparts to 
the proposed two MIR shells are revealed in the PV diagrams. In 
addition, clumps are located along the shells which suggests that 
they are formed within them. In order to estimate the mass and 
kinematic properties of the shell structure, we constructed a cir- 
cular PV diagram which azimuthally averages spectra at a given 
radius from the center. The center was chosen to maximize the 
velocity difference between peaks in the C ls O 7 = 2-1 spectra, 
and is close to the UC Hn regions O and N. Moreover, the circu- 
lar PV diagram (Fig.[3]c) shows a clear elliptical signature which 
implies an expanding motion centered on a common origin. 
Therefore, we conclude that the two MIR shells share a common 
expansion center with an average inner radius of m 2.9 pc (10% 
uncertainty) and an expanding velocity of ~ 5 km s _1 . We can 
reject the p ossibilit y that the s hells are actually molecular rings 
(Bea umont & Williamsl 2009) because two-component spectra 
are seen toward the center of the shells probing their front and 
back sides (Fig. [3]). In the LTE assumption, we can estimate the 
H2 column density of the expanding shells from the 13 CO/C 18 
7 = 2-1 ratio, assuming their optical depth ratio is the same as 
their abundance ratio . We adopted 12 C/ 13 C of 46 and 16 0/ 18 of 
510 dWilson & Roodl [1994: Lang er & Penziasl[T990T) for W49A 
(8.1 kpc to the Galactic center, fawinn. Moran. & ReidL [l992h . 



The excitation temperature of the averaged shell is estimated to 
be about 25 K. The average H2 column density of the shell is 
a 8.5 x 10 22 cm -2 , excluding the contribution from source G 
(i.e., consider the shell density within a radius of 40" in Fig.[3]d). 
Assuming a 10% uncertainty, we calculate an H2 density, n(H2), 
of * 6.7+2.1 x 10 3 cm 3 in a spherical structur e with a filling fac - 
torof 1 and a CO abundance of 2 x 10~ 4 (e.g jLacv el al.. 1994). 
The shell mass (M sh ) can be estimated as 1.9 + 0.7 x 10 4 M©. 
Then the kinetic energy (jM^V^) and the momentum (M s hV s h) 
can be computed to be 4.6 + 3.0 x 10 4S erg and 9.3 + 4.0 x 10 4 
M© km s _1 , respectively, assuming a 20% uncertainty in the ex- 
pansion speed measurement. 



4. Discussion 

4. 1 . Driving mechanisms for the expanding shells 

The reason for the elliptical shape of the two shells is not clear, 
but is probably due to different expansion velocities in the two 
axes, or due t o multiple expanding events. In the interstellar bub- 
ble model of IWeaver et al.l (I1977I) . where the spherical expand- 
ing shells are driven by winds, the age of an averaged expanding 
shell in W49A can be estimated as 3.4 + 0.4 x 10 5 yr from the 
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Fig. 3. Left panel: The linear PV diagrams of 13 CO 7 = 2 - 1 and C ls O 7 = 2 - 1 corresponding to Cut I (E-W) and Cut II (N-S) 
shown in Figure Q] The suggested shell expansion center and source G are marked. The red dashed ellipses indicate the expansion 
speed of 5 km s _1 and the radius corresponding to the two MIR shells marked in FigureQ] The clumps associated with source O and 
source N are marked. (Ia)-(Ib) Black contours represent the main-beam brightness temperatures 7mb from 5% (3 <x for C 18 and 
10 cr for 13 CO) to 35% in steps of 5% of the peak value, and the subsequent contours are running from 40% to 90% in steps of 10% 
of the peak value. (Ila)-(IIb) Black contours represent T MB from 10% (3 <x for C 18 and 10 cr for 13 CO) to 35% in steps of 5% of 
the peak value, and the subsequent contours are running from 40% to 90% in steps of 10% of the peak value, (c) The circular PV 
diagram of C ls O 7 = 2-1. The radius axis is the distance to the shell expansion center, and the temperature is averaged in a given 
radius. Black contours represent the main-beam brightness temperatures of C 18 7 = 2-1 from 0.4 to 2.0 K in steps of 0.2 K. (d) 
The H2 column density of an averaged shell calculated in the LTE assumption. Black contours represent the H2 column densities of 
5.0 x 10 21 cm~ 2 , 1.0 x 10 22 cm -2 , 1.5 x 10 22 cm~ 2 , and 2.0 x 10 22 cnT 2 . The shell with a radius size of 3.3 pc and an expansion 
speed of 5 km s _I is indicated in red ellipses. 
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where E is the mechanical luminosity (E = lM w Vj), and no 
is the ambient H density. We derived an average H2 density of 
~ 2 x 10 3 cm" 3 in W49A by using the 13 CO/C I8 7 = 2 -1 
ratio mentioned above, and we adopted an ambient density «o 
of 4 x 10 3 cm -3 and a stellar wind velocity V K of 2000 km s _1 . 
Hence, a constant mass-loss rate of ~ 1.2 x 10~ 6 Mq yr~' is 
needed to sustain a wind-driven bubble in a size of ~ 2.9 pc, 
which corresponds to the mass loss rate of one O-type star. In 
addition, about 20% of the wind en ergy is channe l ed into the ki- 
netic energy of the swept-up shell dWeaver et all Il977l) . which 
means that a total energy of a; 3.4 + 2.2 x 10 48 erg is provided 
by winds during the expansion. This is consistent with the aver- 
age shell kinetic energy of 4.6 + 3.0 x 10 48 erg. However, given 



the elliptical shape of the two shells, their expansion likely dif- 
fers from that of the interstellar bubble model discussed above. 
Therefore, a reasonable upper-limit of the shell age (f < ^) is 

estimated to be 5.7 + 1.3 x 10 5 yr, assuming either a monoton- 
ically decreasing expansion velocity or a constant velocity. We 
conclude that the ages of these expanding shells are likely about 
3-7 x 10 5 yr. 

On the other hand, a high radiation pressure in Hn regions 
can also provide enough energy to drive an expanding s hell (e.g., 
iKrumholz & Matznerll2009UHomeier & Alvesl (120051) estimate 
a cluster mass of ~ 1 x 10 4 Mq around source O and source 
N by summing 54 stars with masses greater than 20 Mq within 
2.5 pc, which means that this cluster probably hosts ~ 30-50 
Q6 stars with an ionization lum inosity Q of ~ 1 x 10 51 s 1 
(IVacca. Garmanv. & Shulll [l996). This corresponds to a lumi- 
nosity L of 8 x 10 4() erg s~' . The luminosity (L - P^c/ t) needed 
for energy-driven shells can be estimated to be ~ 8 x 10 39 erg 
s _1 , assuming no photon leaking. Therefore, no matter whether 
the expanding shells are driven by winds or by radiation pres- 
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Fig. 4. Upper panels: The left panel shows the large-scale GRS 13 CO 7=1-0 emission of the W49 complex, and the right 
panel shows the 13 CO 7 = 2-1 emission from IRAM 30m in the local area of W49A. The MIR shell structure is indicated by 
black ellipses. The white dashed lines represent two linear PV cuts shown in the lower panels. The white dotted line s represent the 
sugge sted alignments or ejections of molecular gas. Bl ack crosses represent the 850 urn dus t continuum peaks from Matthew s et al] 
(2009). Black triangles denote the UC Hn regions from lDe Pree. Mehringer. & Gossl(fT997h . (III)-(IV) The NW-SE and NE-SW PV 
diagrams corresponding to Cut III and Cut IV in the upper-left panel. The shell with a radius size of 3.3 pc and an expansion speed 
of 5 km s _1 is indicated in red ellipses. The violet curves in a double-cone shape represent the gas ejections in the NW-SE direction, 
and the red dashed line represents the NE-SW direction of the bipolar gas ejections. The contours denote the temperatures from 0.4 
to 4.8 K in steps of 0.4 K. 



sure, just a few massive stars can account for the driving en- 
er gy rather than the whole cluster . A re cent model of W49A 
bv lMurrav. Ouataert. & Thompsonl d2010h demonstrated that the 
dominating radiation pressure from a new born cluster forms a 
bubble which eventually disrupts the whole cloud. In their one- 
dimensional model, the outward force is dominated by the force 
from protostellar jets during the first 3 x 10 5 yr, and radiation 
pressure dominates until the first supernova explodes. However, 
in the early phases, the expanding shell is driven by protostellar 
jets in their model, and it is not clear how those jets can drive 
the shells in W49A. Furthermore, leaking photons may reduce 
the radiation force on the s hell. It is possible t hat m any of the 
embedded sources seen by iHomeier & Alvesl (120051) are stars 
formed after shell fragmentation. We suspect in fact that these 
two shells originate from a common c enter, which is close to th e 
embedded cluster center indicated by Hom eier & Alvesl d2005h . 
and trigger star formation later when expanding. 



4.2. Fragmentation of the shells 



According to Whit worth et al.l (1 19941) . an expanding shell begins 
to fragment at ff rag ment ~ (^Po-M) -1 ' 2 , depending on the ambi- 
ent density po and the shell Mach number At The time scale 
for fragmentation, ffragment, is shorter than the time scale needed 
for the shell to become dominated by self-gravity, f grav ity ~ 
(Gpsh) -1 ^ 2 . The value of fragment is estimated to be » 3.6 x 10 5 yr, 
with a sound speed Ch, of 0.35 km s and a Mach number of 
about 14 inside the shell. The time scale for the shells in W49A 
to fragment into clumps and further collapse to form new stars is 
about a few times 10 5 yr, which roughly corresponds to the time 
scale of massive star formation (Zinne cker & Yorkel feOO?). The 
clumps formed after fragmentation are well separated as shown 
in Figure [3] and the two velocity components at ~ 4 km s -1 and 
~ 12 km s _1 toward the region around source G are perhaps two 
clumps along the line of sight. 
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4.3. Large-scale gas ejections 



Table 1. The projected kinetic energy estimates of gas ejections 



A recent paper by iMatthews et alj ((2009) reveals the distri- 
bution of dust clumps in the large-scale environment of the 
W49 complex, including o n e UC Hn region G43.24-0.05 
dKurtz. Churchwell. & Woodl 1994b and one hydroxy ! maser 
source OH43. 13+0.04 (iBecker. White. & Procton.ll992h . Many 
dust clumps are aligned toward the region near source O and 
source N, which are close to the center of our expanding 
shells (Fig. |4j). Furthermore, the PV diagrams de rived using the 
Galactic Ring Survey (GRS: lJackson et al.L l2006) 13 CO 7=1-0 
emission show an "outflow structure" which may be related to 
the expanding shells. For example, the PV diagram of Cut IV 
shown in Figure [4] indicates a bipolar outflow between G43.24- 
0.05 and OH 43.13+0.04 with the center of the expanding shell 
structure located roughly symmetrically between the two clumps 
seen in mm emission. Other evidence is the detection of UC 
Hn regions (Fig. |4| toward the shell expansion center as well 
as the fact that the filaments seen in 13 CO 7 = 2-1 have a 
radial geometry relative to the shell expansion center (although 
those filaments may be also related to W49A North or source 
G). Still further evidence of gas ejections comes from a large- 
scale gas distribution in the NW-SE direction suggesting that 
widely-distributed molecular gas and UC Hn regions are present 
together with the double-cone structure shown in Figure [4] (III). 

Without taking the inclination into account, we can esti- 
mate the projected kinetic energy (as a lower-limit) of the NE- 
SW gas ejections by assuming that they occurred at the same 
time when the shells started expanding (i.e., 5.7 x 10 5 yr as 
an upper-Umit) and b y adopting the clump mass estimates of 
Matt hews et alJ d2009h . In the northeastern ejection, the clumps 
G43.24-0.05 and G43.224-0.038 have a total projected ki- 
netic energy of 1.3 + 0.6 x 10 5l) erg; in the southwestern ejec- 
tion, W49A Southwest, OH43. 13+0.04, submillimeter sources 
G43. 108+0.044 and G43.099-0.050 have a total projected ki- 
netic energy of 1.1 ± 0.3 x 10 50 erg, assuming a 20% uncer- 
tainty in the mass measurements. The projected kinetic energy 
estimates are summarized in Table 1. The projected kinetic en- 
ergy for both directions is roughly equal, and we infer a total 
projected kinetic energy of 2.4 ± 0.7 x 10 50 erg. This result de- 
pends on our assumption about the age of the outflows. If the 
gas ejection time scale were seven times longer than the age of 
expanding shells, the kinetic energy for the shells and the out- 
flows would be similar. In either case, the total energy released 
by shell expansion and gas ejection has a magnitude of 10 50 erg. 
Since we only estimated the kinetic energy for the NE-SW gas 
ejection, the total kinetic energy of the gas ejection event in the 
W49 complex may be larger. It is worth noting that G43.24- 
0.05 i s located at a similar distance (1 1.7^^ kpc: I Watson et all 
120031) as W49A, and probably hosts a newly born star produc- 
ing the UC Hn region. In addition, the nearby source G43.224- 
0.038 is likely related to G43.24-0.05. The nature of the OH 
maser source OH43. 13+0.04 is not clear, since it is not asso- 
ciated wi th any of the Infrared A s tron omical Satellite (IRAS) 
sources dBecker. White. & ProctorLll992h . Additionally, the na- 
ture of the submm sources G43. 108+0.044 and G43.009+0.050 
are also unknown. 

Th e magnetic field measurement by iBrogan & Trolandl 
(1200 11) shows a gradient in the line-of-sight component of the 
magnetic field, Z?i os , in the NE-SW direction crossing the expan- 
sion center, where a very weak B\ os is derived from the 7 km s _I 
Hi component. For a radially increasing magnetic field, a weak 
Bios is expected in the shell expansion center (Vlsr ~ 7-8 km 
s~'). The fact that Z?i os increases in the northeastern and south- 
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G43.224-0.038 


12.6 


2.6 x 10 3 


21.7 ± 


4.9 


1.2 + 0.6 


G43.24-0.05 


15.6 


1.6 x 10 4 


26.8 ± 


6.0 


11.4 + 5.6 


W49A Southwest 


6.2 


3.5 x 10 4 


10.7 ± 


2.4 


4.0 ± 2.0 


OH43. 13+0.04 


12.2 


8.0 x 10 3 


21.0 ± 


4.7 


3.5 ± 1.7 


G43. 108+0.044 


15.9 


1.9 x 10 3 


27.4 ± 


6.1 


1.4 + 0.7 


G43.009+0.050 


18.0 


2.0 x 10 3 


31.0 ± 


6.9 


1.9 + 0.9 



a The distance to the shell expansion center without a projection cor- 
rection. 

b The gas mass of clumps adopted from Matth ew's et alj (|2009). 

c The ejection velocity of clumps assuming a constant ejection ve- 
locity in 5.7 x 10 5 yr. 



western directions may be related to the bipolar outflow in the 
NE-SW direction. In addition, another strong B\ os peak in the 
northwestern direction at 4 km s _1 could be related to the NW- 
SE gas ejection with a double-cone shape. It is not clear what 
kind of event can cause large-scale gas ejections with a kinetic 
energy of a few times 10 50 erg. A supernova explosion is un- 
likely since such an explosion would destroy the whole cloud 
with a kinetic energy of ~ 10 51 erg. A stellar merger of two 30 
M(T) stars can releas e a ki netic energy as high as ~ 2 x 10 50 erg 
(iBallv & Zinnecketi [2005). However, there is no evidence for 
such an event, and it may be impossible to find conclusive evi- 
dence. Thus, understanding the driving source requires examin- 
ing the region around the shell expansion center in detail. 



5. Conclusions 

The morphology of the molecular line emission and the asso- 
ciated velocity signatures are consistent with the shell struc- 
ture seen in the GLIMPSE images. The expanding shells in 
W49A are suggestive of triggered massive star formation in just 
~ 10 5 yr. There is evidence for dense clumps coincident with 
the shells, and the two velocity components toward the region 
around source G may be two distinct clumps along the line of 
sight. The expanding shells are likely to have a common origin 
close to source O with an expansion speed of ~ 5 km s _1 and 
an age of 3-7 x 10 5 years. However, it is not clear what is the 
driving mechanism of the shell expansion. The average mass of 
the expanding shells is ~ 2 x 10 4 Mq with a kinetic energy of 
~ 10 49 erg, which is probably powered by a few massive stars 
in the center of the cluster, and the expanding shells are either 
driven by strong winds or by radiation pressure. In addition, evi- 
dence of the gas ejections is present in the large-scale GRS 13 CO 
7=1-0 and dust continuum data. Both ejections appear to have 
the same origin as the expanding shells with a total energy of a 
few times 10 50 erg. The properties of the central source driving 
the expanding shells and gas ejections need further investigation. 
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